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Abstract 

We study the gravitational collapse of Type I matter field in N dimensional spacetime 
with radial pressure p r as a function of r. We find that for a given smooth initial data set 
satisfying physical requirements, naked singularities exist for spacetime dimensions N = 4 and 
5 while for N > 6 Cosmic Censorship Conjecture withholds its ground. We, also, study the 
collapse with linear equation of state and find that, similar to dust collapse with appropriate 
choice of initial data naked singularities occur in all dimensions. 
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I. Introduction 

In the recent paper P R. Goswami and P. S. Joshi pointed out that the occurrence of naked 
singularity (NS) as end state of dust collapse beginning from smooth initial data with phys- 
ically relevant restrictions, can be removed when one goes to higher dimensional space-time. 
This has restored Cosmic Censorship Conjecture (CCC) (articulated by Roger Penrose [2]) in 
dust collapse for sufficiently high dimension of the space-time. Many other cases of the dust 
collapse have been studied by various authors 0111 . 

A lot of importance is being given to the physically relevant models such as Type I matter 
fields that include most of the known physical forms of matter like dust, perfect fluid, etc. The 
authors in [SJEIEI have discussed the role of initial data in spherically symmetric gravitational 
collapse for Type I matter fields. Recently, in the paper |Hj A. Mahajan et . al have shown 
the restoration of CCC in Einstein cluster spacetime for spacetime dimensions N > 6 with 
non-zero tangential pressure pg and radial pressure p r = 0. 

F.I. Cooperstock and et al. jH] have shown that if the radial pressure p r is positive , r > 
then non-central singularity is covered by the horizon irrespective of the sign of tangential 
pressure pe ■ 

In recent years, special attention is given to higher dimensional gravitational collapse 
since the development in string theory indicate that gravity may be a higher dimensional 
interaction. Hence, it would be interesting to know the status of CCC in higher dimension 
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spacetime. These studies may help in possible appropriate mathematical formulation of the 
conjecture or assist in providing any possible proof of CCC. 

We take up the problem of Type I matter field collapse with radial pressure p r as a function 
of r and arbitrary tangential pressure pg . The apparent horizon for the evolving collapse has 
been studied and it is found that the final state of collapsing cloud results in the black hole 
(BH) for the dimension of the spacetime N > 6. 

Introduction of equation of state has importance as it makes the collapsing models more 
physically relevant. Docs the linear equation of state p = kp stimulates the outcome of NS/BH 
phases as end state of gravitational collapse? We study, this case in section V. 

In section II, we briefly summarize the field equations and state the conditions on the 
initial data under which the collapse evolves. The evolving field equations with p r = p r (r) 
have been illustrated in section III. In section IV, the behaviour of the apparent horizon is 
studied. The case with linear equation of state has been analyzed in section V. Finally, we 
make some concluding remarks. 



II. Field equations and initial data 

The general spherically symmetric metric describing N dimensional space-time geometry 
within the collapsing cloud can be described in comoving coordinates {t, r, 9, 4>) by 



ds 2 = _ e 2u{t,r) dt 2 + e mt,r) dr 2 + R 2 ^ r)dn 
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is the metric on an (N — 2) sphere. The stress-energy tensor for Type I field in a diagonal 
form is given by ^H] 

T\ = -p, T r r = p,., T e e =T+=p e (2) 

The quantities p, p r and pg are the energy density, radial and tangential pressures respectively. 
We take the matter field to satisfy weak energy condition i.e. the energy density measured by 
any local observer be non-negative, so for any vector V z , we must have, TikV l V k > which 
means p>0 ; p + p r >0 and p + pg>0 . 

Einstein field equations for the metric (1) are described by 
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where we have defined G(t, r) = e 2 ^(R') 2 and H(t, 
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The arbitrary function F(t, r) represents the total mass in a shell of collapsing cloud of 
comoving radius r. The weak energy conditions imply F > 0. The regularity at the initial 
epoch t = ti is preserved by F(ti,0) — i.e. the mass function should vanish at the centre of 
the cloud. The space-time density singularity is due to R — and R' = 0, the later one is due 
to shell-crossings and can possibly be removed from the space-time hence we consider 
here only a physically relevant singularity R = known as shell-focusing singularity where all 
matter shells collapse to a zero physical radius. The scaling freedom available for the radial 
co-ordinate r is being used to introduce the function v(t,r) by the relation 

R(t,r)=rv{t,r), (8) 

this relation is due to the definition v(t,r) = R/r [T2|. We have v(t i: r) = 1; v(t s (r),r) = 
and for collapse v < 0. The time t = t 8 (r) corresponds to the shell-focusing singularity R = 0. 
The six arbitrary functions of the shell radius r as given by v(ti, r) = v (f), ip{ti, r) = i/j {t), 
R(U,r) = r, p(ti,r) = p (r), p r (U,r) =p r „(r), pe{U,r) = p 6o {r) evolve the dynamics of the 
initial data prescribed at the initial epoch t = i,-. We have a total of five equations with seven 
unknowns, namely p,p r ,pg,v,ip, R and F giving us the freedom of choice of two functions. 
Selection of these two free functions, subject to the given initial data and the weak energy 
condition above, determines the matter distribution and the metric of the spacetime and 
thus, leads to a particular time evolution of the initial data. The existence and uniqueness of 
solution of the system of field equations with above mentioned initial data has been discussed 
by Joshi and Dwivedi The solution continues to exist in the neighbourhood of the 

singularity given by R = 0. 



III. Radial pressure p r = p r (r) 

We consider radial pressure P r as a function of r. The two allowed free functions p r and 
v(t,r) are chosen as follows: 

p r =p r {r), v{t 7 r)=c{t)+ V {R) (9) 

Now, let us see, how the Einstein's field equations react to this choice which will decide the 
end state of collapse. As p r is a function of r alone, integrating equation we obtain 

nt.r) = - (N 2 _ P [* N _ 2) + *(r) 

where z(r) is another arbitrary function of r. We need to choose z(r), 

z(r) = [4/(iV — 1)(N — 2)]r N ^ 1 p r (r) so that the mass function takes the form 

2 p r (2r»-i - R"-i) 

(JV-l)(JV-2) 1 j 

that satisfies the regularity conditions at the initial epoch F(ti, r) = 2 p r r N ~ 1 / {N — 1)(JV — 2) 
and F(ti,0) = 0. From equation l|10(l . it is evident that non-negativity of radial pressure 
maintains non-negativity of mass function throughout the gravitational collapse ie. F > 
provided p r > 0. Next, equation © becomes 

_ (N - 1) p r (2r N - 2 - R N - 2 R') + p' r (2r JV ~ 1 - R*- 1 ) 
P ~ {N-1)R N ~ 2 R' ( ' 

At the initial epoch t = U, 

Po=Pr + j^— [ rp r >0 (12) 
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, an all important relation between energy density, radial pressure and pressure gradient. In 
the limit of approach to the singularity ie. as r — > and t — ► t g 

9= ^f=i , (13) 

the density blows up at the singularity. Equation ( ) can be written in the form 

G = e 2 r >^d(r) (14) 

with d(r) = 1 + r 2 b(r) where b(r) is the energy distribution function for the collapsing shells. 
Using equation QlOj l in equation (JSj, we obtain 

(P + Pr) v,r R'R=(N-2)(p e - Pr )R' -Rp' r (15) 

The tangential pressure blows up in the limit of approach to the singularity as the density 
diverges at the singularity and we observe that all of initial data is not independent at the 
initial epoch t = ij 

(p + Pr) r),Rr = (N-2)(pe-Pr)-rp' r . (16) 

Using the initial data prescribed at the initial epoch t = ti to evolve the collapse, we can 
integrate equation (|16(l and obtain, 

V{R) = / ^7 — x ~dR (17) 

Jo R{Po+PrJ 

so, velocity distribution function ?y(i?) is determined by smooth functions p r and pg. Using 
equations I|1U|) and l|14[l . equation (jfj) can be put in the form 



R (N-2),2 R = a{t) e V + r 2 b{r)]R N-3 e 2 v _ R N-S + 2 ^T^_^1 {l8) 

where a(t) is a function of time and we choose a(t) = 1 by suitable scaling of time co-ordinate. 
We define 

e 2')(A) _ 1 

h(R)= R2 = 2g(R) + 0(R 2 ) (19) 
Using this definition, equation l|18|l becomes 



v {N-i)/2 i= _ &v L N -Z h r r)e 2 v + V N-I h f rv ) + JPr V N % ) (20 ) 

y w (N — 1)(N -2) 

Integrating above equation, we obtain 

t(v,r)= / (21) 

Jv e ^v N - s e 2 %(r) + v N ^h(rv) + 2 -^^gLJ2 

where the variable r is treated as a constant. Further, the time taken for the central shell to 
reach the singularity is given by 

ts = / 22 
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ancU So is well defined provided [v N - 3 b(0)+v N - 1 h(0) + [2 p r (0) (2-w Ar - 1 )]/[(iV- l)(iV-2)]] > 
0. The Taylor expansion of t(v, r) around the center r — is given by 



t(v,r) =t(v,0)+r Xp (v) + O{r 2 ) 

The time taken for other shells close to center r — to reach the singularity as t 
determined from 

t s (r) = t So + r Xp {0) + 0{r 2 ) 

and we have 



(23) 
t s can be 

(24) 
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It is clear that Xp( v ) depends on the functions p r (0),h(0) and 6(0) and these functions, in 
turn, depend on the given set of density, tangential pressure and velocity function v. For any 
constant v surface, we have dv = which implies v Xp( v ) = — v ' ■ 
Now, equation 1)20(1 takes the form 



fo(«)] ^(M 



where 



ZJ(0,») 



,JV-3 



6(0) 



.,N- 



ih(0) 



2p r (0) (2-^-i) 



(26) 



(27) 



(iV- l)(JV-2) 

From equations l|16() and (|12() . it is evident that all of initial data is not independent. Hence, 
we conclude that b(r), radial and tangential pressures p r and pe will be sufficient to form 
initial data set for the collapse. Now, we assume that the initial data functions 6(r), pg(r,v) 
and p r (r) are at least C 2 functions and that satisfy the regularity conditions at the center 
r = ie. p' ro (0) = p' e (0) = and p To (0) - p 6o (0) = 0. 
Using equation (|26fl . we can write 



Xp (0) = hm o x^) = lim^-^ 



(28) 



Since the function D(0,v) > and v E [0, 1] and if the above limit exists, then the sign of 
X P (0) as positive or negative shall depend on the behaviour of v' . Hence, v' > may give 
rise to Xp(0) > as limit of sequence of positive terms can not be negative. If Xp(0) = then 
we need to consider the higher derivative of t(v,r) at r — for similar analysis. 



IV. Study of apparent horizon 

The end state of gravitational collapse, in terms of either a naked singularity or a black 
hole, is determined by the causal structure of non-spacelike curves in the neighbourhood of 
singularity The singularity is visible if there exist future directed families of non-spacelike 
curves emanating from the singularity that reach the far away observers in the future, and 
which have past end point at the singularity. Otherwise event horizon forms early enough 
to cover the singularity and thus, the end state of collapse is a black hole. To determine 
this characteristic of singularity, we study the behaviour of apparent horizon which is the 
boundary of trapped surfaces as the collapse evolves. This boundary of the trapped region of 
the space-time within the collapsing cloud is given by the equation, 
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and which is the equation for the apparent horizon, for details, we refer the reader to 
Firstly to know about the NS /BH phases for any r > along the singularity curve t = t s (r) . 
For example, from above equation along the singularity curve (which corresponds to R = 0), 
for any r > 0, the mass function F(t, r) takes the form 



F(t,r)= Cl (2 



,N- 



for a specific choice of r, whereas the physical area radius R{— C2v) — > as v — > where Ci 
and C2 are constants. Therefore, trapped surfaces form that cover the neighbourhood of the 
center before the development of singularity. 

Using equation 1)29(1 ■ the apparent horizon v^ah is obtained through the equation 



^-r + (iV-l)(iV-V-3_ 2 = , 



2r 2 p r 



(30) 



We know from above that all the points r > on the singularity curve are already covered 
and the behaviour of apparent horizon is being studied close to the central singularity at 
r = 0, R = 0, therefore, the values of r in this analysis must be sufficiently small. With this 
background, the equation (|30|l is a polynomial equation in v. For dimensions N = 4 and 
N = 5 the VNah are determined as positive roots of above polynomial equation and we obtain 
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For iV > 6, the process of determining the solution of equation (|30|l goes tedious. Here, 
we argue that by choosing r sufficiently small, v N_1 in equation (|30|1 can be neglected . For 
example, for N = 6, r = 0.01 and p r (r) = 1 + r 2 + r 3 + r 4 + r 5 + r 6 , equation (|30|l 
gives u = 0.02714509001 . With this, we observe that for N = 6 the term /~ 3 [(N - 
1)(N - 2)}/[2r 2 p r ] = 1.999999985 which is very close to 2, while the contribution from the 
term v N ~ x = 0.1473860672 x 10~ 7 is negligible, hence, we can do way with the later term. 
Therefore, for N > 6, we have 



VNah 



4r 2 P r 



(N - l)(N -2) 



(33) 



Now, we discuss the nature of central singularity r — 0, t — t s (0) through the behaviour of 
causal paths emerging from the singularity and reaching out to a far away observer in future 
or otherwise. Using equation <|21|l . the equation of apparent horizon in (t, r) is written as 
follows 

t ah {r)=t s {r)~B(r) (34) 



where 



B(r) 



v (N-3)/2 dv 
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,N-3 
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Since the behaviour of apparent horizon is being studied close to the central singularity at 
r = 0, R = 0, the upper limit of integration in the above equation is small. Hence, the 
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integrand in above equation can be expanded in a power series in v and keeping only leading 
term, we obtain 

t ah (r) = t So + r Xp (0) - Jj^j (Pr)^ (vNah)^ (36) 

For N — 4, 5 and if Xp(0) is non-zero positive, then in above equation the second term 
dominates over the last negative term so that the apparent horizon curve is increasing as we 
move away from the singularity, therefore the singularity is naked. On the other hand for 
N > 6, the negative term in equation H36|) starts dominating and thus, advancing the time 
of trapped surface formation. Hence for N > 6, we have a black hole solution. These results 
coincide with the results obtained in [8] with radial pressure p r = 0. 

These results can be verified by plotting curves for x P (0) = 0.034, P r = 1 + r 2 + r 3 + r 4 + 
r 5 + r 6 on small range of r say to 0.075 and tah(r) can be plotted with range of 0.9994 to 
1.0024. Further, with another value of X P (0) = 0.5, improving the range of r and tah(r), we 
observe that curves for N = 4, 5 escape to infinity. 



V. Collapse with linear equation of state 

In this section, we aim at the question, whether the choice of equation of state p — kp, 
k 6 [0,1] contributes through k in the development of BH/NS phases. R. Goswami and P. 
S. Joshi have studied the case of an isentropic perfect fluid with linear equation of state in 
four dimensional spacetime, wherein, they pointed that the occurrence of NS/BH evolving 
from regular initial data depend on the choice of rest of the free function available [T^. The 
stress-energy tensor for Type I field in a diagonal form is given by |10| 

T\ = -p, T r =T\ = T\=p (37) 

The quantities p and p = p r = pg are the energy density and pressure respectively. We take 
the matter field to satisfy weak energy condition which implies p>0; p + p>0. Linear equation 
of state for perfect fluid is 

p(t, r) = k p(t, r) where k e [0, 1] (38) 

Einstein field equations for the metric (1) are derived as 

(N-2)T _ 1 {N-2)P 
' 2R N - 2 R> k 2R N ~ 2 R 

^ = -^ I pn(p)]' (40) 
R'G - 2RG is' = (41) 

G ~ H = 1 -^- ( 42 ) 

As before, we have R(t, r) — r v(t, r) and the choice of physically reasonable initial data 
is maintained herein as in section II. For further details, we refer the reader to [Sj- A general 
mass function for the cloud can be considered as 

F(t,r) = r N ~ 1 M(r,v) (43) 

where M{r, v) is regular and continuously twice differentiable. Using equation Q43|) in equations 
lj39"|) . we obtain 

_ jV-2 (N — 1)M+ r[M. r + M„ v'] _ (N — 2)M„ 
P ~^2~ X v N ~ 2 (v + rv') ~ 2k v N ~ 2 ( ' 
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Then as v — > , p — * oo and p — + oo i.e. both the density and pressure blow up at the 
singularity. At the initial epoch, the regular density distribution takes the form 

N —2 

Po{r) = — 2— [rM(r, 1), r + (N - l)M(r, 1)]. (45) 

From equation l|44|). it is clear that p = p(r,v) and therefore, v' = f(r,v). We rearrange 
equation (|4"4"|l as follows 

(N - l)kM+ krM t T + Z(r, u)M „ = (46) 

where 

Z(r,v) = (fc + l)rV + v (47) 
Equation 146|l has general solution of the form, 

S(X,Y) = (48) 
where X(r,v, M) and Y(r,v, M) are the solutions of the system of equations |14j . 

dM dr dv 



(N-l)k kr Z 



(49) 



Equation (|48|l has many classes of solutions but only those classes of solutions should be 
considered which satisfy energy conditions, which are regular and gives p — » oo as v — > 0. 
This means the energy conditions and equation of state p = kp isolate the class of functions 
M(r,v) so that the mass function !F(t,r), the metric function v(t,r) and the function b(r) ( 
to follow) evolve as the collapse begins according to the Type I field equations. 
On integrating equation 140(1 . we obtain the general metric function, 

p(r,v)=-^-j[ln(p)] (50) 

Define a suitably differentiable function A(r, v), A(r, v ) v = v' /R' so that equation (|5L)|) takes 
the form 

<51) 

At the initial epoch, we have 

^■^'-wM^ ,52) 

whereas using equation 1441) . the relation between the function Mand A, at all the epochs is 
given by 



A(r,v) tV R' = - 



(fc + 1) 



ln( _(iV-2)M, 



2k v 



N-2 



(53) 



From above, it is clear that M(r, v) is a decreasing function with respect to v and if we consider 
a smooth initial profile for the density in such way that density gradient vanishes at the center 
then we have A(r,v) = rg (r,v) where g {r,v) is another suitably differentiable function. 
Also, the use of definition of A(r, v) in equation (|4*T)> yields 

G(t,r) =d(r)e 2rA (54) 

where d(r) is another arbitrary continuously differentiable function of r. We can write 

d(r) = 1 + r 2 b(r). (55) 
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where b is the energy distribution function for the collapsing shells. Define a function h(r, v) 

as 

'ZrA i 

h M) = ^— (56) 
and substituting this equation, together with l|55|l and (|54|l in equation 1421) . we get 



v (N-3)/2 i = _ p -k/(k+l) ^JV-3 h ( r) „) + b „iV-3 e 2rA + M (57) 

where negative sign is chosen since, for the collapse, v < 0. Integrating the above equation , 
we have 

t(v,r) = / 58 

J v p -k/(k+l) v / w A r -3 h ( r . j + bu Ar-3 e 2rA + M 

In above equation, the variable r is treated as a constant. Expanding t(i>, r) around the center, 
we get 

t(v,r) = t(v,0) +r X (v) + 0{r 2 ) (59) 

where the function 

dt _ 1 f 1 v( N -V/ 2 BA0,v) 



X{V) dr 



2j v B(0, vf/ 2 



-dv. (60) 



and 

B(r, v) = /9 - fe /( fe + 1 ) y/v N - 3 h(r, v) + bv N - 3 e 2rA + M(r, v) (61) 

-1 v (N-3)/2 dv 

-mw- (62) 

t s (r) = t(0, r) = t So + rx(0) + 0(r 2 ) (63) 

Now, it is clearly seen that the value of x(0) depends on the free functions b(0), .M(0, v) 
and h(0,w), which in turn, depend on the initial data at the initial surface t — U. Thus, a 
tangent to the singularity curve t = t So is completely determined by the given set of density, 
pressure, velocity function v and function b. Further, from equation H57J) . we can write 

v {N - 3)/2 v' = X (v)B{0,v) + O(r). (64) 

Hence, 

x(0) = lim x(v) = Hm . (65) 

«->o v^o B[Q,v) 

The existence of the limit depends on the cumulative effect of all the terms present in above 
equation but since B(0,v) > as v — > provided M(0,0) > and v S [0,1], therefore 
positive or negative sign of x(0) absolutely depends on v'. So, if v' > then %(0) > and 
otherwise. If x(0) — then we will have to take into account next higher order non-zero term 
in the singularity curve equation and do a similar analysis. Thus, we reckon that value of 
x(0) depends on the initial data profiles of density and pressure (through the free function 
M(r, 1)), the energy function b(r) and the other free function h(r, 1) (for determination of 
the metric function ip). Secondly, its sign depends on the nature of gradient of v. 

Now, we assume that these free functions to be at least C 2 and satisfy other physical 
requirements. Let us choose the function M (r, v) as follows 

M(r, v) = 1 + (rv) + (rv) 2 + (rv) 3 + 

Now, the equation of apparent horizon r n-& = 1 give rise to a polynomial equation in v 

v N - 3 -rv N - 2 =r 2 . (66) 
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The apparent horizon V ^ah for space dimensions N = 4, 5 and 6 are derived as follows 



^ 1 ± Vl - 4r 3 
2r 

next for iV = 5, we have 



A = -108r 4 + 8 + 12a/3 \/27r 4 - 4r 2 



which is an imaginary term for sufficiently small values of r, < r < y/ 4/27. Hence, apparent 
horizon may occur at some value of r > ^/ A/27 for a chosen initial data of M(r,v). So, for 
such value of r, we can write 

A 1/3 2/3 1/3 



> rA 1 / 3 



r 



and 



1 V3S , 1 J -185T + 6(12) 1 /3 r 5T 2 / 3 + 24(12) 2 / 3 r 4 S'- 18VSTV3 



4r 12r 12r V ST 1 / 3 



where T = r 2 (9 + \/-768r 5 + 81) and 5 = ^ 3Tl/3 + (212)V^rT 2 / 3 + 8(12) 2 / 3 r 4 ^ 

The nature of central singularity r = 0,t = t s (0) as NS/BH is studied through the behaviour of 
geodesies emanating from the singularity to reach a asymptotic observer in future or otherwise. 
Using equation l|58[) . the equation of apparent horizon in (t, r) is written as 

t ah (r)=t s (r)-B 1 (r) (67) 

where 

rVah v (N-3)/2 d 

Bi(r) = / (68) 

Jo p -fc/(fc+l) ^/„Af-3 h ( rj + hv N-3 e 2rA + M 

The upper limit of integration in the above equation is small as discussed earlier. Using the 
value of p from equation (|44|l and expanding the terms in the integrand in a power series in 
v and keeping only leading term, we obtain 

Mr) _ t . (r) _ '<~-"" r r't::^ m 



(1 - rw)-( 1 - 2 ™)/ 2 

where we have put m = k/(k + 1) and < m < 1/2. Further, since r is sufficiently small, 
therefore, we can do away with the terms 0(rv) 2 , hence for N > 4, we obtain 



[(JV — 1)(AT — 2)1 T 
t ah r = tso + rxO + ~ 7^ - 



r(l - 2m)Viv a/l 



(V/Va/i 



1 -(iV+l-2m)/2 

iV+l-2m AT — 2m — 1 

(70) 

The second term in above equation dominates over the last for sufficiently small values of r 
and for any value of k (the parameter of equation of state) and independent of any choice 
of dimension of the spacetime. Hence, causal paths emerging from the singularity reach out 
to an asymptotic observer at infinity, thus naked singularities appear in all dimensions for 
sufficiently small values of r. The results for dust case ie. with m = agree with the results 
obtained in 
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VI. Conclusion 



We have considered the radial pressure p r as a function of r, due to which the mass function 
F(t, r) becomes the function of both r and v, while the tangential pressure remains arbitrary 
through out the analysis. Through the study of apparent horizon, it is found that NS forms at 
the center of the collapsing cloud for dimensions N = 4 and N = 5, whereas for TV > 6 there 
exists a critical value of x(0) below which a BH solution persists. This statement supports 
the understanding that gravity may be a higher dimensional interaction. This generalizes the 
results obtained in 

Next, we have analyzed the effects of the equation of state on the formation of appar- 
ent horizon in continual gravitational collapse, the parameter k involves in determination of 
the curve t — t a h{r) and assist in strengthening the out going causal path to reach out to 
an asymptotic observer at infinity. It is analyzed that naked singularities exist in all the 
dimensions, k = is the usual case of dust collapse. 
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